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ABSTRACT 


Atmospheric  density  fluctuations  induce  phase 
perturbations  that  degrade  the  spatial  coherence  of  a  laser 
beam  propagating  through  the  atmosphere.  These 
degradations  spread  the  laser  beam  and  alter  the  centroid 
and  intensity  profile  stochastically.  Turbulent  conditions 
arising  from  various  environmental  situations  and 
meteorological  phenomena  are  found  at  virtually  all  levels 
of  the  atmosphere.  A  substantial  fraction  of  the  optical 
turbulence  along  a  vertical  path  arises  from  the  surface 
heat  flux  within  the  first  100-200  meters  above  the  ground. 

This  thesis  seeks  to  measure  and  analyze  these 
turbulent  layers.  A  high  frequency  acoustic  echosounder 
was  developed  to  analyze  atmospheric  turbulence  within  the 
first  100-200  meters  above  the  ground  with  extremely  high 
accuracy.  The  echosounder  design  incorporated  a  25 
element,  square,  planar  array  housed  within  a  lead  and  foam 
insulated,  acoustic  shroud.  This  acoustic  echosounder  was 
employed  to  collect  real-time,  low  level  atmospheric  data 
at  two  sites.  This  data  should  prove  useful  in  determining 
whether  the  performance  of  laser  and  electro-optical 
systems  could  be  increased  by  raising  sensors  and 
transmitters  well  above  the  local  turbulent  layer  through 
the  use  of  towers  or  other  such  structures. 
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A  coherent  laser  beam  propagating  through  the 
atmosphere  is  very  susceptible  to  numerous  turbulent- 
dependent  processes.  As  electromagnetic  waves  transit  the 
turbulent  regions,  atmospheric  irregularities  randomize  the 
amplitude  and  phase  of  the  wave.  A  high  resolution 
acoustic  profiler  was  designed  and  developed  to  quantify 
the  altitude  dependence  of  these  atmospheric  disturbances. 
Acoustic  energy  is  more  susceptible  to  these  irregularities 
than  electromagnetic  waves.  Historically  acoustic  sounders 
have  proven  to  be  the  best  means  of  providing  high 
resolution  atmospheric  turbulence  profiles  within  the  first 
500  to  1000  meters  above  the  ground.  Acoustic  echosounders 
are  frequently  used  to  detect  and  measure  atmospheric 
density  and  velocity  fluctuations  resulting  from  wind 
shears,  convection,  and  temperature  inversions. 

Presently  two  optical  parameters,  the  spatial  coherence 
length  ( ro )  and  the  isoplanatic  angle  ( 80 )  are  measures  of 
the  perturbation  of  electromagnetic  waves  propagating 
through  the  atmosphere,  and  are  measured  very  reliably  by 
optical  systems  developed  by  Walters  [Refs.  1  and  2]. 
Although  these  systems  measure  a  line  integral  of  the 
atmospheric  optical  turbulence  with  high  accuracy,  a  major 
drawback  of  each  system  is  that  no  provision  is  provided 


for  determing  the  height  of  the  atmospheric  disruptions. 
If  these  disturbances  are  found  to  exist  very  near  the 
surface,  it  may  be  possible  to  negate  their  effect  by 
elevating  sensors  and  transmitters  or  altering  the  surface 
heat  flux  to  suppress  the  generation  of  turbulence. 

This  thesis  deals  with  the  design,  construction  and 
implementation  of  a  high  frequency,  acoustic  echosounder 
which  will  accurately  analyze  the  atmospheric  density 
fluctuations  within  approximately  the  first  200  meters  of 
the  atmosphere.  As  this  project  is  a  product  of  the 
research  and  efforts  of  two  students,  the  work  was 
appropriately  divided.  My  particular  task  was  concentrated 
in  the  actual  design  and  hardware  development  of  the 
echosounder.  My  colleage,  LT .  Michael  Wroblewski  [Ref.  3], 
dedicated  his  efforts  toward  software  development  and  data 
display,  analysis  and  storage. 

The  echosounder  developed  and  discussed  herein  utilizes 
a  25  speaker,  square,  planar  antenna  array  enclosed  within 
an  acoustically  absorbent  shroud  comprised  of  lead  and  foam 
insulation  inside  a  sturdy  plywood  framework.  This  antenna 
array  was  employed  as  both  the  transmitter  and  receiver  of 
the  acoustic  signal.  A  computer  controlled  function 
generator  transmits  an  acoustic  pulse  through  the  antenna 


array . 


Additional  electronic  components  incorporated 


within  the  overall  echosounder  system  are  used  to  amplify, 
digitize,  display,  analyze  and  ultimately  store  the  return 


wir'WA 


signal.  Over  time  the  received  signal  data  is  used  to 
produce  a  quantitative,  three  dimensional  map  of  the 
atmospheric  turbulence. 

Acoustic  echosounders  have  been  developed  and  in  use 
for  many  years  and  have  proven  to  be  a  valuable  probe  for 
analyzing  the  structure  and  dynamics  of  the  lower 
atmosphere  [Refs.  4  through  7].  Devices  similar  to  ours 
have  been  used  to  obtain  profiles  of  the  atmospheric 
density  and  temperature  fluctuations  [Refs.  8  and  9],  Our 
device  excels  by  utilizing  a  high  speed  HP217  computer  to 
control  and  monitor  the  echosounder.  This  enables  us  to 
obtain  real-time  data  collection  and  analysis  with  the 
ability  to  store  and  reproduce  the  atmospheric  profile 
plots  at  will.  This  information  should  help  quantify  the 
effects  of  low  level  turbulence  on  electro-optical  systems 
performance.  It  may  then  be  possible  to  negate  these 
turbulent  effects  by  simply  mounting  sensors  and 
transmitters  well  above  the  turbulent  layers  or  by  adaptive 
optics . 

This  thesis  discusses  the  development  of  a  high 
frequency,  5KHz,  acoustic  echosounder,  the  criteria 
involved  in  the  actual  designs  of  both  the  acoustic  array 
and  the  acoustically  absorbent  enclosure,  and  the  selection 
and  function  of  all  the  associated  component  parts  of  the 
echosounder  system.  Finally  the  operational  success  of  the 
system  will  be  discussed  and  the  output  data  analyzed. 


9 


II.  BACKGROUND 


Acoustic  echosounders  probe  the  atmosphere  ty 
transmitting  a  pulse  of  acoustic  power  which  i  r- 
subsequently  scattered  back  from  the  atmosphere  by 
temperature  and  velocity  inhomogeneities.  The  echosounder 
( echosonde )  equation,  often  referred  to  as  the  radar 
equation  in  meteorology,  is  used  to  determine  the 
backscattered  acoustic  power.  This  equation  is  summarized 
by  Neff  in  Reference  10  and  is  based  upon  the  work  of 
Tatarski  [Ref.  11]  and  Little  [Ref.  5], 

Pr  =  Er  [  Pt  Et  ][  exp(-2aR)  ][  ao(R.f)  ][  ^ci  ][  AGR- ^  ] 

where 

Pr  is  the  electrical  power  returned  from  a  range  R, 

Pt  is  the  electrical  power  transmitted  at  frequency  f, 
Er  is  the  efficiency  of  conversion  from  acoustic  power 
to  electrical  power  by  the  transducer, 

Et  is  the  efficiency  of  conversion  from  electrical 
power  to  acoustic  power  by  the  transducer, 

exp(-2aR)  is  the  round  trip  power  loss  due  to 
attenuation  where  a  is  the  average  attenuation  (meters  1  > 
to  the  scattering  volume  at  the  range  R  (meters), 

oo(R,f)  is  the  scattering  cross  section  per  unit  volume 
at  a  distance  R  and  frequency  f, 


c  is  the  local  speed  of  sound  (meters/second), 
r  is  the  pulse  length  (seconds), 

A  is  the  aperture  area  of  the  antenna  (meters2  ), 
R  is  the  range  (meters),  and 


f. 


1  is  the  effective-aperture  factor  of  the  antenna. 

Empirically  measuring  or  calculating  the  values  for  all 
"ther  terms,  one  can  use  this  equation  to  determine 
’'j'R.fi,  the  scattering  cross  section  per  unit  volume;  that 
is  the  fraction  of  incident  power  backscattered  per  unit 
iistance  into  a  unit  solid  angle  at  a  frequency  f.  Based 
jpon  experimental  results,  Tatarski  [Ref.  11]  expresses  the 
acoustic  backscatter  cross  section  per  unit  volume, 
co  R  .  f  '  .  by. 


c ft  R, t  i 


0.0039  k1  /  3 


Ct  2 


To2 


wne  r^ 

k  -  2n/\  is  the  incident  acoustic  wavenumber  at 

wave  1 ength  A  , 

To  is  the  local  mean  temperature  in  degrees  Kelvin,  and 

Ct 2  is  the  temperature  structure  parameter. 

'ombining  this  equation  with  the  echosonde  equation,  one 
obtains  a  volume -averaged  measure  of  Ct 2 . 


11  2  1  Pr 

- - To2  k_1/3  - - - R2exp(2aR) 

U.0039  Er  Et  ct  AG  Pt 


1  1 
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•a 


a 


Til 
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Hall  and  Wescott  [Ref.  12]  calculated  a  beam-shape 


compensation  factor  of  0.40  for  a  piston  source  antenna 


with  a  uniformly  illuminated  square  aperture.  This  value 


is  the  same  as  the  effective-aperture  factor,  G  and  can  be 


substituted  into  the  above  equation.  Approximating  the 


aperture  area  of  the  antenna  to  be  equivalent  to  25  times 


the  aperture  area  of  a  single  speaker  having  a  diameter  of 


7.620  centimeters,  we  get  a  value  of  0.1140  square  meters. 


Combining  these  values  with  the  numerical  constants  in  the 


above  equation,  we  can  simplify  the  equation  for  Ct2. 


1  1  Pr 

Ct2  =  11245  - To2  k-i/3 - - R2exp(2aR) 

Er  Et  ct  Pt 


The  efficiency  of  conversion  factors  for  our 


echosounder  were  determined  with  the  use  of  an  IVIE 


Electronics  IE-30A  Audio  Analyzer  Sound  Level  Meter,  a 


HP3314A  Function  Generator,  a  HP3561A  Dynamic  Signal 


Analyzer,  a  reference  speaker,  the  acoustic  array  and  the 


Naval  Postgraduate  School  anechoic  chamber.  The  function 


generator  was  used  to  transmit  a  1.0  volt,  5.0  KHz, 


sinusoidal  signal  through  the  reference  speaker.  The 


acoustic  array  was  employed  as  a  receiver  and  was  aligned 


along  the  centroid  of  the  main  lobe  of  the  reference 


speaker  signal.  With  the  acoustic  array  in  position,  the 


sound  level  meter  was  used  to  determine  the  total  power 


available  for  reception  by  the  acoustic  array;  that  is,  the 
total  power  actually  ensonified  within  the  solid  angle  of 
the  acoustic  array.  Ten  independent  measurements  were  made 
with  the  sound  level  meter  and  then  averaged  to  yield  a 
power  value  of  73.7  dB.  Using  a  conversion  factor  and  the 
array  aperture  area  determined  above,  the  total  acoustical 
power,  Pa ,  available  at  the  acoustic  array  is  given  by  the 
equation  below. 

Pa  ~  107-37  •  (  1  •  10-12  Watts  /  m2  )  •  (  0.11401  m2  ) 

or , 

Pa  =  2.67  •  10-6  Watts 

The  dynamic  signal  analyzer  was  then  used  to  measure  the 
actual  power  received  by  the  acoustic  array  and  converted 
into  electrical  power.  Four  independent  measurements  were 
recorded  and  averaged  to  yield  a  value  of  4.035  rnVrms  .  The 
acoustic  array  impedance  was  previously  measured  and 
determined  to  be  12.1  ±0.4  Q.  These  values  were  used  in 
the  following  equation  to  determine  Pe ,  the  total 
electrical  power  converted  by  the  acoustic  array. 

(  4.035  •  10-3  Vr  ms  )  2 

Pe  - -  =  1.34  •  10-6  Watts 

12.1  Q 


The  efficiency  of  conversion  factor  for 


acoustic  to 


electrical  power  is  the  ratio  of  Pa  to  Pe ,  or 


3 


— H 
| 

i 


a 


-*,5 

<1 


1.34  •  10-6  Watts 


Er  = 


2.67  •  10-6  Watts 


0.503 


The  value  for  the  efficiency  of  conversion  for  electrical 
to  acoustic  power,  Et ,  is  not  readily  calculatable  with  the 
equipment  at  hand.  However,  based  upon  the  design  of  the 
speakers  used,  the  values  for  Er  and  Et  should  be 


approximately  equal. 


Furthermore,  allowing  for 


experimental  error  especially  in  the  reading  of  the  sound 
level  meter,  the  best  accuracy  for  the  values  of  Er  and  Et 
is  0.5  ±0.1.  Substituting  conversion  factor  values  of  0.5 
into  the  above  equation,  we  can  further  simplify  the 
equation  for  Ct2 . 


To2  Pr 

Ct2  =  44980  k-t/3 - Ra  exp(  2aR) 

CT  PT 


The  above  equation  was  incorporated  into  a  computer 
program  that  provides  the  system  control,  data  acquisition 
and  data  processing  techniques  [Ref.  3].  The  reduced  data 
was  then  used  to  measure  the  temperature  structure 
parameter  as  a  function  of  time  of  day  and  altitude  for 
various  sites. 


A.  SPEAKERS 


In  the  design  of  our  echosounder,  it  was  determined 
that  a  rapid  decay  time  was  required  to  obtain  accurate 
short  range  information.  Weight  restrictions  involved  with 
equipment  transportation  plus  high  efficiency  in  the  2KHz 
to  5KHz  range  led  to  our  decision  to  use  piezo  ceramic 
speakers.  The  Motorola  KSN  1005A  speaker  was  selected 
based  on  these  requirements  and  the  specifications  charted 
in  the  Motorola  catalog  [Ref.  13]  and  reproduced  below  in 
Figures  1  through  3. 

Nominal  Power,  Impedance,  and  Distortion  Ratings 

UJ  401— r 
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Typical  Frequency  Response 


INPUT  2  8V 

MICROPHONE  DISTANCE  1/2  Mtltr 
(1  Will  into  8  Ohms) 


FREQUENCY  (Hz) 

500  1000  2000  5000  10000  20000  40000 

Fig.  2.  Speaker  Frequency  Response 

Dimensions:  KSN  1005A,  KSN  1003A 
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WEIGHT:  75  GRAMS 


Fig.  3.  Speaker  Dimensions 
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These  charts  indicate  that  a  maximum  response  for  our 
speakers  occurs  at  a  resonant  frequency  of  5000  Hertz. 
This  frequency  was  used  as  the  baseline  from  which  all  our 
measurements  were  made. 

Using  the  speakers  in  an  anechoic  chamber,  the  average 
e* 1  voltage  decay  time  was  measured  to  be  approximately  900 
usee  (Fig.  4).  This  decay  time  translates  into  a  sound 
propagation  distance  of  just  over  15.0  centimeters  (at  STP) 
from  the  speakers.  Considering  our  requirements,  this 
speaker  was  ideally  suited  to  serve  our  purpose. 


B.  ACOUSTIC  ARRAY 


Our  next  consideration  was  the  echosounder  array 
pattern.  Ideally,  the  acoustic  sources  should  be  placed 
exactly  one  half  wavelength  apart.  At  a  frequency  of  5000 
Hertz,  this  would  require  spacings  of  3.4  centimeters  (at 
STP)  which  is  physically  impossible  for  the  speakers  we 
have  chosen.  The  closest  possible  spacing  was  7.62 
centimeters  between  sources  after  shaving  off  the  flange  of 
the  horn  (Fig.  3). 

From  Kinsler,  et  al .  [Ref.  14:  p.  195],  the  equation 
for  the  directionality  factor  of  a  simple  line  array  is 
derived  as: 

N 

sin  (  -  kd  sin  0  ) 

2 


1 

sin  (  -  kd  sin  8  ) 

2 

where 

k  is  the  wave  number  (2x/i  ), 
d  is  the  distance  between  sources, 

N  is  the  number  of  sources ,  and 
0  is  the  angle  measured  from  a  line  perpendicular  to 
the  array  to  the  direction  of  interest. 

However.  this  equation  assumes  simple  point  sources 
which  does  not  accurately  describe  the  speakers  we  have 
chosen.  It  was  necessary  to  couple  this  equation  to  the 


directionality  factor  for  a  piston  source  which  is  also 
identified  in  Kinsler,  et  al .  [Ref.  14:  p.  108]. 


D(9) 


2  Ji  (  ka  sin  6  ) 


ka  sin  0 


where 

k  is  the  wave  number, 

a  is  the  radius  of  the  piston  source, 

0  is  the  angle  measured  from  a  line  perpendicular  to 
the  array  to  the  direction  of  interest,  and 
Ji  is  a  first,  order  Bessel  function. 

These  two  equations  were  combined  to  produce  the 
equation  of  directionality  for  a  linear  array  of  piston 
sources,  L(N,8),  by  simple  multiplication. 

L(N, 0  )  =  H(N,0)  •  D( 0 ) 

This  equation  was  incorporated 
T Appendix  A]  and  used  to  generate 
varying  numbers  of  piston  sources 


into  a  computer  program 
linear  array  patterns  for 
(Figs.  5  through  9). 
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Fig.  5.  Three  Element  Array  Fig.  6.  Four  Element  Array 


Fig.  7.  Five  Element  Array 


Fig.  8.  Six  Element  Array 


Fig.  9.  Seven  Element  Array 


It  was  concluded  from  these  plots  of  the  linear  array  beam 
patterns  for  varying  numbers  of  piston  sources  that  five 
linear  elements  would  provide  the  best  combination  of 
forward  directionality ,  sidelobe  suppression,  physical 
size,  and  a  relatively  low  cost.  Then,  in  order  to  enhance 
both  array  efficiency  and  symmetry,  we  settled  on  a  five  by 
five  element  array  design  with  7.62  centimeter  spacing 
between  speakers  in  both  the  vertical  and  horizontal 
directions . 

After  we  verified  the  manufacturer’s  polarity 
designation  for  35  speakers  and  inspected  each  of  these 
speakers  for  physical  defects  and  loose  wiring  connections, 
we  obtained  Lissajous  plots  for  each  individual  speaker. 
Figure  10  represents  a  typical  Lissajous  plot  for  one  of 
the  speakers  incorporated  into  the  five  by  five  element 
array.  Based  upon  the  speaker’s  output  to  input  voltage 
ratios  as  illustrated  by  these  Lissajous  plots,  we  were 
able  to  rank  all  of  our  speakers  by  signal  efficiency.  It 
was  based  upon  this  criteria  that  we  selected  the  25  most 
efficient  speakers  for  use  in  the  five  by  five  element 
array,  placing  the  best  speakers  at  the  center  and 
subsequent  ranking  speakers  further  toward  the  sides  and 


corners . 


Fig.  10.  Lissajous  Plot 

The  25  selected  speakers  were  mounted  in  a  five  by  five 
planar  array  on  a  balsa  wood  insulated  bilayered  sheet 
metal  board.  After  wiring  all  the  speakers  in  series,  we 
surrounded  all  the  electrical  connections  and  speaker  backs 
with  two  3.0  centimeter  layers  of  foam  insulation 
sandwiching  a  1.0  millimeter  lead  sheet.  Then  the  entire 
array  mounting  was  enclosed  in  a  44  by  44  by  5  centimeter 
sheet  metal  box.  This  design  was  chosen  to  suppress 


virtually  all  acoustic  energy  propagating  out  the  rear 
hemisphere  of  the  array,  while  shielding  the  array  from  any 
external  electrical  interference  (Figs.  11  and  12). 


C .  ENCLOSURE 


Acoustic  echosounding  has  proven  to  be  an  extremely 
useful  technique  for  probing  and  analyzing  the  lower 


atmosphere.  In  order  to  utilize  the  acoustic  waves 
transmitted  and  later  received  by  this  remote  sensing 
method  efficiently,  it  is  essential  to  have  an  efficient 
antenna  with  highly  directive  beams  and  strongly  suppressed 
sidelobes.  Antenna  design  becomes  increasingly  more 
important  in  a  noisy  environment  where  noise  pollution 
within  the  sidelobes  may  dominate  the  desired  signal  within 
the  main  lobe.  Hall  and  Wescott  [Ref.  12]  showed  that 
sidelobe  suppression  improved  with  higher  frequencies. 
Their  studies  showed  that  the  measured  90  degree  sidelobe 
suppression  ranged  from  38  dB  at  1  KHz  to  50  dB  at  5  KHz. 
Furthermore,  any  significant  improvement  in  sidelobe 
suppression  could  only  be  obtained  by  surrounding  the 
antenna  with  an  acoustic  energy  absorbing  cuff  or  shroud. 

In  an  effort  to  maximize  our  antenna  main  lobe  to 
sidelobe  power  ratio,  we  intend  to  operate  at  5KHz . 
Additionally,  we  have  designed  an  acoustic  energy  absorbing 
enclosure.  Many  designs  were  considered  based  upon 
previous  research  in  the  field  of  echosounding  [Refs.  15 
through  18],  In  addition,  we  obtained  the  actual  acoustic 
beam  patterns  for  our  array  using  a  computer  program 
written  by  LCDR  Butler  [Ref.  19]  which  we  modified  for  our 
purposes.  This  modified  version  of  LCDR  Butler’s  computer 
program  is  in  Appendix  B.  By  rotating  the  array  in  an 
anechoic  chamber,  we  were  able  to  produce  highly  accurate 
polar  plots  of  tne  array  beam  patterns  (Figs.  13  and  14). 
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Fig.  13.  Polar  Plot  of  5  X  5  Array  Beam  Pattern 
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A  comparison  of  the  theoretical  beam  pattern  for  the 
five  element,  piston  source  array  in  Figure  7  with  the 
experimentally  measured  polar  plot  of  the  5  X  5  acoustic 
array  in  Figure  13  and  the  associated  coordinates 
corresponding  to  the  individual  data  points  illustrate 
significant  agreement.  As  predicted  by  our  computer  model 
[Appendix  A] ,  the  actual  acoustic  beam  pattern  is  comprised 
of  nine  distinct  lobes  in  the  forward  hemisphere.  Both  the 
positions  of  these  lobes  and  their  relative  amplitudes  are 
within  2.0  percent  of  their  predicted  values.  The  computer 
modeled  main  lobe  full  width  angle  of  19.5  degrees  very 
closely  predicted  the  experimentally  measured  main  lobe 
full  width  angle  of  19.0  degrees.  Reflected,  scattered, 
diffracted  and  electronic  noise  all  contribute  to  the 
clutter  near  the  origin  of  Figure  13.  This  is  the  primary 
reason  why  the  minima  don’t  actually  go  to  zero  as 
predicted  by  our  model.  Additional  reasons  include 
mechanical  noise  in  the  device  used  to  rotate  the  acoustic 
array  and  vibrational  noise  in  the  array  itself.  The 
apparent  lobe  in  the  upper  left  hand  quadrant  is  an  anomaly 
seen  in  all  polar  plots  herein,  and  corresponds  to  acoustic 
energy  reflected  off  the  door  and  window  of  the  anechoic 
chamber.  Figure  14  represents  the  acoustic  array  mounted 
at  a  45°  aspect.  A  computer  model  for  this  orientation  is 
a  very  complicated  task  and  was  not  performed  at  this  time. 
However,  this  polar  plot  of  the  acoustic  array  beam 


suppression  strongly  support  the  need  for  further  array 
design  studies. 

Since  it  was  our  aim  to  suppress  all  sidelobes  and 
utilize  the  main  lobe,  we  chose  not  to  taper  our  enclosure 
as  most  previous  researchers  had.  Rather  we  designed  the 
enclosure  based  upon  the  dimensions  of  the  array  itself  and 
the  acoustic  beams  it  generated.  The  height  of  the 
enclosure  is  160  cm.  Allowing  for  the  base  foam  and  lead 
insulation,  the  array  housing  and  the  base  framework 
itself,  the  actual  array  piston  sources  are  positioned  143 
cm  from  the  mouth  of  the  enclosure.  The  length  and  width 
of  the  enclosure  are  both  65  cm.  Again  allowing  for 
construction  materials,  the  inner  dimensions  are  both  48 
cm.  Hence,  the  position  of  the  center  speaker  is  24  cm 
from  any  of  the  four  sides.  These  dimensions  produce  the 
unobstructed  19.0  degree  path  into  which  the  main  lobe  is 
transmitted . 

2  •  tan-i  (  24  cm  /  143  cm  )  =  19.050  „  19.00 

Plywood  was  used  for  the  construction  of  the  enclosure 
and  provided  not  only  a  rigid,  inexpensive  framework,  but 
also  provided  an  impedance  mismatch  which  helped  to 
attenuate  external  noise  interference.  Anticipating  all 
kinds  of  weather  conditions  during  data  collection,  the 
plywood  enclosure  was  first  waterproofed  with  four  coats  of 


marine  varnish.  Grooved  joints,  caulking  and  weather 
stripping  were  also  design  considerations. 

A  millimeter  layer  of  lead  can  suppress  an  acoustic 
signal  by  reflection  as  much  as  40  dB  (Figs.  15  and  16). 
About  7.0  centimeters  of  convoluted  foam  can  suppress  a 
signal  another  3  to  4  dB  (Figs.  15  and  17).  The  primary 
purpose  of  the  convoluted  foam  is  to  attenuate  the 
sidelobes  which  repeatedly  reflect  off  the  enclosure 
interior  at  glancing  incident  angles.  Together  the  lead  and 
convoluted  foam  make  an  extremely  efficient  absorbing 
material  for  use  in  our  enclosure. 


START:  4  000  Hz 
X :  3000  Hz 


RANGE:  -19  cfBV  pA*JSED 


BU:  19.097  Hz 
y: -19.29  dBV 


STOP:  8  000  H 


surfaces  of  the  enclosure  with  two  1.0  millimeter  layers  of 


lead  overlapping  at  all  corners.  Strong  aluminum  brackets 


well  as  to  the  enclosure  base  (Figs.  18  through  20) 


Enclosure  Photograph,  Fully  Assembled 


D.  ADDITIONAL  HARDWARE 

In  addition  to  the  acoustic  array  and  enclosure 
discussed  in  the  previous  sections,  there  are  primarily  six 
more  pieces  of  equipment  essential  to  the  function  of  the 
echosounder.  Below  is  a  brief  description  of  these 
additional  components. 

1 .  HP200  Series  Computer 

The  Hewlett-Packard  (HP)  200  Series  Computer 
includes  a  20  megabyte  hard  disk  drive,  a  floppy  disk  drive 
and  an  associated  printer  and  monitor.  The  HP200  Series 
Computer  is  the  central  control  component  for  the  entire 
echosounder  arrangement.  The  computer  used  is  an  HP217 
programmed  in  Basic  3.0  and  equipped  with  an  Infotek  BC203 
Basic  Compiler  and  an  Infotek  FP210  Floating  Point 
Accelerator  to  enhance  the  speed  of  program  execution.  The 
computer  is  used  to  maintain  timing  and  process  control,  to 
define  all  parameters  for  the  HP3314A  Function  Generator  as 
well  as  to  execute  the  trigger  command  that  produces  the 
echosounder  transmitted  acoustic  pulse.  The  computer  also 
processes  the  data  from  the  acoustic  array  via  the  pre¬ 
amplifier,  bandpass  filter  and  analog  to  digital  converter. 
The  computer  then  conducts  data  reduction  routines  to 
produce  and  display  a  high  resolution  atmospheric  profile. 

2 .  HP3314A  Function  Generator 

The  Hewlett-Packard  (HP)  3314A  Function  Generator 
is  a  multimode  function  generator  capable  of  providing 


sine,  square  and  triangular  wave  functions  as  well  as  any 
desired  waveform  ranging  in  frequency  from  0.001  Hertz  to 
19.999  Megahertz.  The  HP3314A  Function  Generator  was  used 
to  supply  a  pulse  of  an  integer  number  of  sinusoidal  cycles 
of  constant  amplitude  to  the  QSC  Model  1700  Audio 
Amplifier.  Typically  a  pulse  of  100  sinusoidal  cycles  was 
used.  A  constant  frequency  setting  of  5000  Hertz  was  used 
for  all  data  runs. 

3 .  QSC  Model  1700  Audio  Amplifier 

The  QSC  Model  1700  Audio  Amplifier  is  a  high  power 
amplifier  which  can  supply  350  watts  over  an  8  ohm  load. 
This  amplifier  was  used  to  boost  the  output  voltage  of  the 
function  generator  by  a  factor  of  20  from  1.5  volts  to  30 
volts.  The  power  supplied  to  each  speaker  in  the  array 
during  operation  was  then  37 . 1  watts  as  indicated  in  the 
equation  below. 

(  Voltagerms  )2  (  30  )2 

Power  - - -  - =  37.1  Watts 

Impedance  2  •  12.1 

4 •  Pre -amplifier 

The  pre-amplifier  was  designed  and  constructed  by 
Walters  [Ref.  20]  to  supply  a  gain  of  1000  to  the  returned 
signal  and  isolate  the  data  aquisition  system  from  the 
transmitted  pulse.  The  pre-amplifier  may  be  thought  of  as 
a  safety  and  switching  mechanism  for  the  system.  An  LT 


1007  and  an  LT  1037  Operational  Amplifier  were  selected  for 


Bandpass  Filter  was  the  best  filter  available.  It 
performed  adequately  despite  the  large  1  KHz  bandpass. 

6 .  Infotek  AD200  Analog  to  Digital  Converter 

An  Infotek  AD200  12  bit  Analog  to  Digital  Converter 
was  used  to  sample  the  amplified  signal  voltage  and 
digitize  it  for  the  computer  system.  A  typical  sampling 
frequency  used  was  20000  Hertz. 


The  arrangement  of  the  above  equipment  i 
in  Figure  23. 
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Fig.  23.  Echosounder  Equipment  Arrangement 
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E .  SOFTWARE 


The  software  and  the  HP217  Computer  were  responsible 
for  controlling  and  monitoring  every  operational  aspect  of 
the  system.  LT  Wroblewski  [Ref.  3]  provides  an  in-depth 
discussion  of  the  echosounder  software  which  is  also 
summarized  in  Appendix  C. 
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A .  ECHOSOUNDER  PERFORMANCE 

An  analysis  of  the  acoustic  echosounder  output  was 
conducted  to  determine  the  validity  of  previously 
determined  echosounder  parameters.  The  e~ 1  decay  time 
constant  for  the  applied  voltage  of  a  single  speaker  was 
previously  calculated  to  be  900  us  which  corresponds  to  a 
distance  of  15.3  cm.  The  e- 1  decay  time  constant  for  the 
applied  voltage  of  the  entire  acoustic  array  was  determined 
by  the  same  experimental  procedure  to  be  1060  us  which 
corresponds  to  a  distance  of  18.0  cm.  Although  the  25 
speakers  in  the  acoustic  array  transmitting  in  unison 
should  have  approximately  the  same  e~  i  decey  time  constant 
as  the  single  speaker,  it  is  not  surprising  that  the  array 
e- 1  decay  time  constant  is  larger.  Using  ten  time 
intervals,  an  expectation  of  99.99  percent  decay  of  the 
transmitted  signal  should  occur  within  a  distance  of  180.0 
cm  or  just  under  2.0  meters.  Reviewing  typical  echosounder 
traces  and  Ct 2  plots  [Appendix  D] ,  it  is  evident  that  such 
a  short  recovery  time  is  actually  quite  accurate.  The 
recovery  distance  is  typically  found  to  be  on  the  order  of 
9  meters  or  53  ms.  This  blind  zone''  is  not  uncommon  in 
most  radar  systems  and  is  in  part  a  return  signal  from 
ground  clutter  surrounding  the  acoustic  array  and  its 


enclosure  immediately  following  pulse  transmission.  An 
analog  to  digital  intensity  factor  of  4000  is  assigned  to 
the  heavily  shaded  region  of  the  recovery  zone.  Since  the 
computer  software  has  generated  this  value  as  the  product 
of  the  actual  analog  to  digital  "count"  value  and  the  range 
from  the  acoustic  echosounder,  this  value  of  4000  must  be 
divided  by  the  range  at  the  recovery  zone  boundary  to  give 
us  an  actual  analog  to  digital  "count"  value.  This  result 
is  then  multiplied  by  an  experimentally  measured  conversion 
factor  of  3.34  •  10-®  voltsr ms /count .  This  result  may  then 
be  used  in  the  determination  of  the  total  electrical  power 
received  and  converted  by  the  acoustic  array.  The  value 
for  the  total  electrical  power  received  by  the  acoustic 
array  is  calculated  in  the  equations  below. 

Voltagerms  =  (3.34  •  10-®v/count)  ( 4000m* counts )  (9m)-i 
Voltagerms  -  1.48  •  10~5  volts 

Power  =  ( Voltager ms  ) 2  •  (Array  Impedance)-1 
Power  =  (1.48  •  10-®  volts)2  •  (12.1  Q)-1 
Power  =  1.82  •  10-ii  Watts 

Using  this  received  electrical  power  value  and  the  known 
transmitted  acoustic  power  value  of  37.1  Watts,  we  can 


solve  for  the  e- 1  decay  time  constant  of  the  applied 


voltage.  The  product  of  the  total  acoustic  power  received 
and  the  efficiency  of  conversion  factor  (  Er  =0.50  ), 
yields  the  total  electrical  power  received.  Hence,  the 
total  acoustic  power  received  is  3.64  •  10" 11  Watts. 
Although  the  distance  of  9  meters  corresponds  to  a  time 
interval  of  53  ms,  this  value  is  actually  inflated  by  the 
system  pulse  length  of  20  ms.  Compensating  for  this,  the 
actual  recovery  zone  time  interval  is  33  ms.  The  e~  2  decay 
time  constant  of  the  power  is  equal  to  the  e- 1  decay  time 
constant  of  the  applied  voltage.  These  values  are 
calculated  below. 


3.64  •  10-11  Watts 


37 . 1  Watts 


e-2t/r  or, 


-  2  •  33  ms 

x  - - -  2.56  ms  , 

In  (6.37  •  10-12) 

This  e-i  decay  time  constant  of  the  applied  voltage 
corresponds  to  a  distance  of  43.5  cm  and  is  about  twice  the 
previously  determined  value. 

Consideration  was  also  given  to  the  construction  of  a 
hexagonal  planar  array  rather  than  a  square  planar  array 
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based  upon  the  astounding  sidelobe  suppression  for  the  five 
by  five  element,  square  array  at  a  45°  orientation  as 
illustrated  in  Figure  14.  To  investigate  the  potential  of 
such  a  close  packed  design,  two  small  seven  element 
hexagonal  arrays  were  constructed.  The  first  array  design 
utilized  a  close  packed  spacing  of  7.62  cm,  while  the 
second  array  design  employed  a  3/2  wavelength  or  10.2  cm 
spacing.  Polar  plots  of  the  beam  patterns  for  both  these 
hexagonal  array  designs  were  made  by  rotating  the  arrays  in 
an  anechoic  chamber  with  the  main  diagonal  both  horizontal 
and  vertical.  In  both  cases  the  first  sidelobe  suppression 
for  both  close  packed  orientations  was  down  by  nearly  15 
dB,  while  the  first  sidelobe  suppression  for  both  3/2 
wavelength  spacing  orientations  was  only  down  by  about  7 
dB.  Furthermore,  the  ratios  of  the  size  of  the  main  lobe 
to  those  of  the  sidelobes  is  much  greater  for  the  close 
packed  design  than  they  are  for  the  3/2  wavelength  spacing 
arrangement  for  both  the  vertical  and  horizontal 
orientations.  The  polar  plots  of  the  array  beam  patterns 
for  both  close  packed  and  3/2  wavelength  spacing  hexagonal 
designs  are  illustrated  in  Figures  24  through  27.  These 
figures  strongly  warrant  further  study  of  the  hexagonal 
close  packed  array  design. 
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Polar  Plot  of  Close  Packed  Hexagonal  Array 
(Horizontal  Diagonal) 
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Polar  Plot,  of  3/2  Wavelength  Hexagonal  Array 
(Horizontal  Diagonal) 
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B.  SITE  EVALUATION 


Echosounder  data  was  collected  at  two  different 
locations.  The  primary  data  collection  site  was  the  upper 
roof  of  Spanagel  Hall  at  the  Naval  Postgraduate  School, 
Monterey,  California.  This  site  was  chosen  simply  for 
convenience  of  data  collection.  Data  gathered  at  this 
location  is  believed  to  represent  the  California  coast 
during  the  spring  near  sea  level.  The  second  data 
collection  site  chosen  was  in  the  vicinity  of  the  24  inch 
telescope  at  Lick  Observatory,  San  Jose,  California.  This 
site  is  located  atop  Mt .  Hamilton  at  an  altitude  of 
approximately  5700  feet  and  nearly  20  miles  inland  from  the 
coast . 

These  two  data  collection  sites  represent  areas  of 
differing  atmospheric  air  pressures,  water  vapor  pressures, 
local  temperature  ranges,  and  local  wind  velocity  ranges. 
These  characteristics  all  play  important  roles  in  effecting 
the  local  atmospheric  turbulent  conditions  and  thereby  the 
atmospheric  structure  parameter,  Ct2. 

In  addition  to  collecting  echosounder  data  at  Lick 
Observatory,  simultaneous  measurements  of  the  isoplanatic 
angle  (9o)  and  spatial  coherence  length  ( ro )  were  made  with 
systems  developed  by  Walters  [Refs.  1  and  2],  A  basic 
knowledge  of  these  two  systems  is  necessary  to  understand 
the  correlation  procedures  made.  The  isoplanatic  angle 


( 9o )  is  primarily  an  upper  atmospheric  measurement  which 
indicates  atmospheric  disruptions  at  a  range  of  two  to  15 
kilometers  above  the  Earth’s  surface.  On  the  otherhand, 
the  spatial  coherence  length  (ro)  is  a  uniform  measure  of 
the  effects  of  the  entire  atmospheric  blanket  on  coherent 
light  transmission.  A  close  comparison  of  the  acoustic 
echosounder  data,  the  isoplanatic  angle  (80)  measurements 
and  the  spatial  coherence  length  ( ro )  measurements  should 
give  us  an  indication  of  the  turbulent  contributions  from 
the  lower  and  upper  troposphere  as  well  as  the  stratosphere 
above . 

Based  upon  isoplanatic  angle  ( 80 )  and  spatial  coherence 
length  ( ro )  measurements  at  Mt .  Wilson  near  Los  Angeles, 
California,  a  strong  correlation  between  the  isoplanatic 
angle  ( 80 )  and  spatial  coherence  length  ( ro )  measurements 
occurs  if  the  low  altitude  boundary  layer  contribution  is 
sufficiently  small.  This  strong  correlation  helps  to 
reinforce  the  overall  description  of  the  atmosphere  at  the 
time  of  data  collection.  Figures  28  and  29  graphically 
illustrate  the  atmospheric  measurements  made  at  Mt.  Wilson 
on  2  April  1987  by  Walters.  The  strong  correlation  between 
the  isoplanatic  angle  ( 0o  )  and  the  spatial  coherence  length 
( ro )  is  especially  evident  between  the  hours  of  0700  and 
1300  universal  time.  The  close  tracking  of  these  two 
measurements  during  this  time  interval  indicate  that  the 


upper  atmospheric  conditions,  as  measured  by  th- 
isopianatic  angle  (do),  are  dominating  the  entire 
atmospheric  profile  as  measured  by  the  spatial  coherence 
length  i ro ' .  Unfortunately,  the  existence  or  non-existence 


if  any  turbulent  surface  effects  cannot  be  asf 


ied  by 


the  employment  of  the  above  two  systems  alone.  However, 
employment  of  these  two  systems  together  with  the  acoustic 
echosounder  should  enable  us  to  produce  a  complete 
atmospheric  profile  with  strong  correlation  between  all 
three  atmospheric  measurements. 

On  9  and  10  April  1987,  all  three  systems  were  operated 
at  Lick  Observatory.  Again  good  correlation  between  the 
i.-epLanatic  angie  (Ho  )  and  the  spatial  coherence  length 
i  ro  )  measurements  was  noted.  In  addition,  strong 
correlation  between  the  spatial  coherence  length  i ro  ) 
measurements  and  the  acoustic  echosounder  data  was 
conclusively  present.  A  comparison  of  the  isopianatic 
angie  (Hu  )  and  the  spatial  coherence  length  fro  i 
measurements  in  Figures  30  and  31  shows  good  correlation 
between  the  two  atmospheric  parameters  especially  during 


the  0830  to  1130  time  interval  on  9  Apri L  1987.  However, 
during  subsequent  hours  the  isopianatic  angle  <  h(,  i 
measurements  remain  high  'approximately  12  urad)  indicating 


1  a  1 1  v 


y  calm  turbulent  conditions  in  the  upp-w' 


mosphere .  while  the  spatial  coherence  length  <  ro  i  vili 


I 


drop  sharply  after  1100  universal  time  indicating  dominant 
and  increasing  lower  atmospheric  turbulence.  This  increase 
in  the  low  level  turbulence  should  be  evident  in  the 
echosounder  data  commencing  around  1200  universal  time 
(0400  local  standard  time)  on  9  April  1987.  A  comparison 
of  the  echosounder  data  in  Figures  32  through  34 
illustrates  this  increase  in  the  local  surface  turbulence. 
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Fig.  30.  Isoplanatic  Angle  Measurements,  Lick  Observatory 
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The  strong  correlation  between  the  echosounder  data  around 
0400  to  0500  local  standard  time  and  the  spatial  coherence 
Length  ( ro )  measurements  at  1200  universal  time,  combined 
with  the  lack  of  correlation  between  the  isoplanatic  angle 
(0o )  and  the  spatial  coherence  length  ( ro )  measurements 
indicate  that  the  lower  atmospheric  and  surface  turbulent 
layers  are  dominating  the  atmospheric  profile  during  this 
time  period. 

Data  collected  on  10  April  1987,  again  illustrate  the 
strong  correlation  between  the  three  atmospheric 
measurements  made.  A  comparison  of  the  isoplanatic  angle 
(0o  )  and  the  spatial  coherence  length  ( ro )  measurements 
'.luring  the  time  interval  of  U700  and  1300  universal  time 
indicate  steady  turbulent  conditions  in  the  upper 
atmosphere  and  greatly  varying  turbulent  conditions  at 
L'.'wer  atmospheric  levels.  This  is  evident  in  Figures  85 
and  36  by  the  consistent  values  of  the  isoplanatic  angle 
i 0o  i  during  the  time  period  compared  with  the  steady 
increase  and  eventual  decline  of  the  spatial  coherence 
Length  ( ro  )  values  during  the  same  time  interval.  The 
trace  variations  in  Figure  36  during  the  hours  of  0800  and 
1300  universal  time  are  indicative  of  a  period  of 
diminishing  lower  atmospheric  or  surface  turbulence 
foi Lowed  by  the  onset  of  an  increasingly  turbulent  period 
nr.und  1200.  This  turbulent  trend  is  distincti vely 
.  uppcrf^d  by  the  echosounder  data  in  Figures  37  through  44. 
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Fig.  44.  Echosounder  Trace,  Lick  Observatory 


I'he  Large  200  mm  spat ia  l  coherence  lengths  around  1  1  "u 
universal  f  ime  are  consist  ent  with  the  low  tu  rbu  1  once 
evident  in  the  echo  sounder  profiles  around  0  3o0  local, 
s  *•  arid  a  rd  t  i  me  . 

drily  acoustic  echosounder  measurements  were  made  on  the 
upper  roof  of  Spanagei  Hall  at  the  Naval  Postgraduate 
'd-hooi.  Data  runs  on  26  and  27  April  1987,  include  bo*:h 
echosounder  data  and  the  associated  Ct 2  plots.  Of 

particular  interest  is  the  plot  of  11:30  on  the  27th  of 
April  in  Figure  45,  which  shows  the  maritime  boundary 
inversion  layer  around  1U0  m  being  perturbed  by  convective 
piumes  at  Lower  altitudes.  Sixteen  additional  echosounder 
output.  traces  were  included  in  Appendix  D  to  exhibit 

typical  atmospheric  activity.  Many  of  the  plots,  such  as 
the  14:45.  15:30,  and  the  17:30  of  26  April  and  the  loom, 

11:15,  13:00,  14:15.  15:15,  16:15,  and  the  17:  lb  of  2 . 

dpri L,  have  convective  plumes  which  are  prevalent  whenever 
a  heat  flux  between  the  surface  and  atmosphere  exists.  The 
plots  of  18:15  and  19:15  on  the  26th  of  April  and  of  13:00, 
18' 45,  and  20:00  on  the  27th  of  April  clearly  show  the 

passage  of  the  neutral  event  which  is  encountered  when  the 

at  m*~.  spheric  and  surface  temperature  difference  becomes 
negligible.  Finally,  the  plot  of  16 : 30  on  26  April  can  be 
as  S'  -c  i  a  t.--d  with  strong  winds  which  exhibit  a  somewhat 
uniform  return  for  all  altitudes  across  the  entire  !‘ 
m  i  n.u  re  in  he  rva  I  of  the  trace  . 


V .  CONCLUSIONS  MD_BECQMMENDAXIQNS 


A  high  frequency,  5kHz,  phased  array  echosounder 
constructed  to  measure  low  level  turbulence  appears  to  work 
within  the  10  to  135  meter  altitude  range.  When  using  the 
device,  detailed  profiles  of  the  short  range  atmospheric 
density  fluctuations  were  obtained.  The  profiles  were 
found  to  correlate  very  well  with  the  measurements  of  the 
isoplanatic  angle  ( 00  )  and  the  spatial  coherence  length 
I ro  )  during  periods  of  simultaneous  operation.  This  short 
range  echosounder,  when  used  in  conjunction  with  the  other 
atmospheric  measuring  devices,  is  an  invaluable  tool.  The 
graphical  output  provides  us  with  a  more  complete 
description  of  the  atmosphere  that  can  be  used  to  calculate 
On 2  ,  the  atmospheric  index  of  refraction  structure 
parameter 

Further  software  developments  such  as  the  incorporation 
of  Doppler  and  Fast  Fourier  Transform  routines  into  the 
echosounder  program  will  determine  the  radial  velocity 
profile  of  the  return  signal.  Such  an  improvement  will 
determine  the  velocity  of  the  probed  air  masses  and  plot 
these  echosounder  traces  as  a  function  of  color  intensity. 

Other  areas  of  further  reasearch  include  the  design  and 
testing  of  different  array  patterns.  One  such  pattern,  the 
hexagonal  array,  is  already  undergoing  tests.  As  the 


speaker  array  evolves,  it  is  also  evident  that  the 
enclosure  must  follow  suit  to  accomodate  the  new  array 
design  and  associated  beam  patterns.  Further  improvements 
in  software  routines  are  also  inevitable.  It  is  always 
desirable  to  store  on  a  disc  all  data  collected  at  a  site. 
Presently,  the  floppy  disc  capacity  is  inadequate  for  the 
storage  of  more  than  two  hours  of  data.  Furthermore,  if 
FFT  and  Doppler  routines  are  to  be  added,  the  computational 
speed  of  the  computer  will  inhibit  the  pulse  repetition 
rate.  A  special  data  acquisition  technique  called  "Direct 
Memory  Access"  may  then  have  to  be  introduced  into  the  HP 
computer  to  allow  simultaneous  data  collection  and 
processing. 

Finally,  this  product  may  be  used  in  various 
applications  not  already  addressed.  One  such  use  may  be  to 
measure  windshear  at  airports.  Other  researchers  have 
expressed  interest  in  using  the  echosounder  to  measure 
arctic  atmospheric  conditions  during  meteorological 
surveys.  Assuredly,  as  the  device  evolves  into  a  more 
compact  and  highly  efficient  instrument,  its  range  of 
application  will  continue  to  expand. 


APPENDIX  A 


LINEAR  ARRAY  PROGRAM 

! ANTENNA  BEAM  PATTERN  FOR  AN  N  ELEMENT  LINEAR  ARRAY 

GCLEAR 

R1-.0381 

Fr=5000 

Ka*( Ri*2*PI*Fr  1/343 

INPUT  "ENTER  THE  #  OF  DRIUERS  “  ,N 

0- .0752 

Kd=( 2*PI*0»Fr )/343 
Kd2=Kd/2 . 

Num_decades=4 

Sea  1 s* 1 0“Nun_de cades 

Thetastep=.001 

GCLEAR 

UIEWPQRT  20,80,20,30 
WINDOW  -20,105,-2,1 10 

AXES  1 00/ <  2*Num_decades ) , 1 00/ ( 2 . *Num_decades ) ,0 ,0 ,20 ,20 
MOUE  10,103 

LABEL  “PATTERN  FOR  “iN;“  ELEMENTS" 

CSIZE  3 
FOR  1-1  TO  5 

MOUE  -15,<2S*< 1-1  >  >-1 
M-< 1 0-<  <  S-I  )*10  )  ) 

LABEL  M 
NEXT  I 
LDIR  PI/2 
MOUE  -15,40 
LABEL  "DECIBELS" 

LDIR  0 

FOR  Theta-Thetastep  TO  PI/2  STEP  Thetastep 
Arg=Ka*SIN<  Theta  ) 

Arg2=2*FNJ 1 ( Arg  ) 

Arg2=Arg2/Arg 
Kds=Kd2*SIN( Theta  ) 

Htheta=SIN(Kds*N )/< N*SIN(Kds >  ) 

Htheta=Htheta*Arg2 
Htheta=Scale*Htheta*Htheta 
Htheta=LGT(  Htheta )/Num_decades 
X-Htheta*SIN< Theta  ) 

Y-Htheta»CQS( Theta  ) 

PLOT  X  *  1 00 , Y* 1 00 
•  PRINT  Theta  .Htheta 
NEXT  Theta 


453  OEF  FNJ1 (X  ) 

4S0  !  Jl  BESSEL  FUNCTION  10  MAY  ,1984;  DLU 

470  IF  X>3 .  THEN  GOTO  Xgt3 
480  X3=X/3 . 

490  X3=X3*X3 

500  Fx=.5+X3*(-.56249985+X3*( . 2 1 093573+X3*< - . 03954289+X3* 

( .004433! 9+X3*( -.00031761 +X3* .0000! 109  >  )  ) )  ) 

510  J 1 =Fx*X 
520  RETURN  Jl 
530  Xgt 3 :  X3=3./X 

540  Fx=.7978845S+X3*( . 00000 1 SS+X3* ( .0! 659667+X3*<  . 000 1 7 1 0S+X3* 

(-.002495! 1 +X3*<  .001 13G53-X3*. 00020033  ) )  )  )  ) 

550  Tx=X-2.35Sl9449+X3»( . 1 2499G 1 2+X3*( . 00005G5+X3*( - . 00G37879+ 

X3«( . 00074348  +  X3*< .00079824-X3*.0002916G  >  > )  )  ) 

5G0  J 1 =Fx*COS( Tx  )/SQR( X  > 

570  RETURN  Jl 
580  FNEND 


APPENDIX  B 


LCDR  BUTLER’S  BEAM  PATTERN  PROGRAM 


10  PRINT  "**  BEAM  PATTERN  COMPUTATION  AND  PLOT  PROGRAM 
20  PRINT  "»*i************»**"*****-*-***********»-*»-**««*»« 
30  ! 

40  1 

50  PRINT  "«****  SET_UP  ****•" 

60  ! 

70  ! 

30  PRINT  CHR$  (10) 

90  PRINT  "THIS  PROGRAM  CONTROLS  THE  FOLLOWING  APPARATU 
100  PRINT  “TO  MEASURE.  DISPLAY.  PLOT  AND  STORE  THE" 

110  PRINT  "DIRECTIONAL  CHARACTERISTICS  OF  A  SINGLET." 
120  PRINT  "DOUBLET.  TRIPLET.  AND  QUADRUPLET  ARRAY.” 

130  ! 

140  PRINT  "  INSTRUMENT  ADDRESS  " 

150  PRINT  "  _  _  " 

160  PRINT  "  HP_3478  MULTIMETER  (VOLTS)  7 24 
170  PRINT  "  HP_3478  MULTIMETER  (RES)  73? 

180  PRINT  "  DATA  AQIJ  IS  I T  ION/CONTROL  709 

190  ? 

200  DISP  "INITIALIZE  INSTRUMENTS" 

210  ! 

220  REM  ***  PRESET  AND  INITIALIZE  INSTRUMENTS--* 

230  OUTPUT  709  ; ”RS"  !  RESET 
240  OUTPUT  723  ;"F3R5T1"  »  FREE  RIJN  30K0HMS 
250  OUTPUT  724  :"F2RAT1"  !  FREE  RUN  AIJTORANGE  A i;  VU!  * 
260  ! 

270  DISP  "INITIALIZE  PROGRAM- 
230  ! 

290  OPTION  BASE  1 

300  COM  ANGLE  (150).  CALC_ANGLE  < 1 SD  >  .  NAMI  5 '  :'i 
310  DEG 

320  PRINTER  IS  701 
330  INTEGER  I.C.L.N 

340  DIM  R<  360)  .VOlT  <  360  )  .OHMS  Da  :m(  ’s;  »  ,  <  <  • 

350  MVOLT=0 

360  MASS  STORAGE  IS  " . DA  T  A 1 " 

370  DISP  "NUMBER  OF  DATA  POINT,  DESIRE:  " 

380  BEEP  10.300 
390  INPUT  M 

400  DISP  "ENTER  DISC  DATA  FILE  NAMt  " 

410  INPUT  NAMES 

420  DISP  "WANT  NEW  DATA  (Y/NW 


n 


4  30 
440 
450 
4b  0 
4  70 
48(1 

4  4 1  j 

50  0 

5  1  0 
500 


BEEP  10.300 
INPUT  Q$ 

IE  QS="Y"  THEN  GOTO  480 
IF  Qt ="N"  THEN  GOTO  300 
GOTO  4  30 

REM  *»*  DEEAIJL  T  PARAMETERS  *** 

MU'. T  INSERT  FOR  LINES  380,  400.  4 !  ij  AND  ylo 
PRIOR  TO  START  OF  RUN 

IREATE  NAME'S.  360.  16 


uj  ( 

.40  ASSIGN#  1  TO  NAMES 


550 

NAME  5 

^ "Rijn  i" 

5b  0 

F  UL  i.  S 

CALF -son  •  •_ 

0(>  IN  MVA< 

5  7  0 

LONG  - 

360/10000  * 

'  UNVER'  OHMS  TO 

DFURFF  5 

58'J 

SF  '  UP  :  •  SE  T 

UP  MEASUREMENT 

f ARAMF  TF  K',  •** 

’  30 

ft?  I  N  T 

X 

XI 
& 

F.IJO 

PR  IN  T 

"'.graph  i: 

5_t  •  NAME  S 

b  T  0 

RF  M  • 

•*  MOD  IF  i  IN 

RRf  HEEDING  L  INF 

UHF  N  NE  LF  GGAKY 

*-> , 1 ' ) 

D  I  SP 

"MANOAl  ,  1  FNTFR  DEGIRFD  F  RF  N 

FOR  SPf  AK  F  R  DR  I 

k  ;o 

DISP 

•*?  N  TR  0  U 

NO  f  HAN(.F  " 

r,4l) 

DI  .P 

’  F  N  T  F  R  '  ’  i. 

C  HANOI  TIC'!" 

a  50 
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920  DISP  "WHAT  IS  DATA  RUN  DESIRED'?" 

930  BEEP  10.500 
940  INPUT  NAMES 
950  ASSIGN#  1  TO  NAMES 
960  FOR  1*1  TO  N 

970  READ#  1  ;  0HMS_DATA< I ) ,VOLT ( I ) 

9.90  NEXT  I 

990  ASSIGN#  1  TO  * 

1000  GOTO  860 
1010  1 
1020  * 

1030  PLQT_ROUTINE:  1  ***  PLOTS  POLAR  DATA 
1040  ! 

1050  * 

1060  PLOTTER  IS  705 
1070  GRAPH 
1080  GCLEAR 
1090  PEN  1 

1100  LIMIT  60 .210 .20 . 180 


AND  ’  •■Ml  PR  I  HR  in  MAR"  'Jf  RON 

i«R  l.  u  <•  -U  mH  i'll'  •!  -  ’J" 


■X-? 


1370 

1380 

1390 

MOO 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1430 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1610 

1615 

1620 


MOVE  L .  0 

FOR  1=0  TO  360  STEP  5 
DRAW  L*COS  < I) ,L*SIN  (I) 

NEXT  I 
NEXT  L 
LINE  TYPE  1 
PEN  UP 
PEN  2 

FOR  1=1  TO  N  !  BEGINS  BEAM  PATTERN  COMPUTATION 
IF  VOLT ( I)>MVOLT  THEN  MVQLT=VOLT( I ) 

NEXT  I 

FOR  1=1  TO  N 

K  =  OHMS_DATA< I) /I  0000*360 
K1=K 

IF  OHMS_DATA< I)<10  THEN  K=K1 +360/295 
IF  OHMS_DATA< I >  >1 0500  THEN  GOTO  1610 
VECTOR=20*LGT  <VOLT< I)/MVOLT)+40 
LINE  TYPE  1 

IF  VECTOR< 1  THEN  VECT0R=O 
PRINT  I, K. VECTOR 
X ( I ) =VECTOR*COS  <K ) 

Y< I ) =VECTOR*SIN  <K) 

PLOT  X < I ) , Y  <  I ) 

NEXT  I 
PEN  UP 
RETURN 
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APPENDIX  C 

ECHOSOUNDER  SOFTWARE 

The  driving  force  of  this  system  is  found  within  the 
HP200  Series  Computer.  The  software  and  computer  system 
were  responsible  for  controlling  and  monitoring  every 
operational  phase  of  the  hardware  components.  The  program 
which  accomplishes  this  task  was  entitled  'ACRDR  and  is 
listed  in  Reference  3.  This  program  was  modeled  after  a 
program  written  by  Walters  [Ref.  20],  but  each  program 
performs  a  distinctly  different  computational  task.  This 
program  was  written  in  HP  Basic  4.0  and  was  compiled  by  an 
Infotek  BC20  i  Basic  Compiler  to  enhance  the  spe^d  t 
execut ion . 

The  program  ACRDR  is  easily  broken  down  into  a  number 
it  blocks  and  subroutines  which  perform  spe.-if: 

■  perations.  These  sections  are  outlined  in  a  rl  w.-har* 
i  Pie  4b).  Each  block  is  straightforward  in  its  purr 
a n d  the  program  is  designed  to  be  as  helpful  to  the  user  a- 
possible.  As  a  prologue  to  the  actual  program  o.-,de  ,  ’  ii-r- 

l  -  a  listing  of  all  the  program  variables  with  t  r.  >  ’ 
i-snption  of  their  use  Du  -h  a  listing  tamii’.  ar:  •*  • 
is-r  with  t  tie  imput  at  i  ns  ♦  •  ■  be  made  a  n  1  .a  i  f  r  .•  :  '• 
lui  ’ k  r-t  »•  r°n'>5  f  ■  r  any  f  ut  ure  m  d  l  f  i  •  •  a »  :  n ’ 


s?NI 


The  initialization  procedures  which  set  all  the 
parameters  necessary  for  data  collection  follow  the 
variable  definition  section.  These  parameters  control  such 
things  as  the  contrast  between  background  and  return  signal 
and  the  setting  of  the  computer’s  internal  clock.  All 
options  of  the  HP3314A  Function  Generator  are  decided  by 
the  user  at  this  phase.  Another  initialization  involves 
preparing  the  Infotek  AD200  Analog  to  Digital  Converter  for 
operation.  Finally,  internal  arrays  are  dimensioned, 
integer  variables  are  defined  and  the  computer  function 
keys  are  redefined  to  suit  the  "ACRDR"  program. 


Initialization 
Procedures  (or: 

HP  331 4A 

AD  2C0  A/D  Convener 
Internal  Clock 
Constants 


Computation  ot 
Parameters 
Speed  of  Sound 
Offset 
Range 

Pulse  Length 
Attenuation  Coefficient 


P'ecare  Sc'een  !cr 
Ca:a  D'So!ay 


T rgger  Pulse  from 

HP  331  4A 


Send  Cj  Plot  to  the 
Printer 


Compute  Values  of  CT“' 
and  Plot 


There  are  certain  options  built  into  the  program  to 
allow  a  user  to  change  various  aspects  of  operation.  The 
function  keys  allow  the  user  to  change  the  local 
temperature,  sample  frequency  or  intensity  factor  during 
program  execution.  Additional ly ,  the  user  can  quit  or 
restart  the  program,  print  the  partial  trace  on  the  screen 
or  elect  to  save  a  trace  on  a  floppy  disc.  The  save 
routine  is  invoked  for  the  subsequent  15  minute  interval 
after  the  appropriate  function  key  is  depress,  d.  Saving  a 
future  trace  may  be  ackward,  especially  if  a  -ser  would 
like  t o  keep  an  interesting  trace  which  is  presently  on  the 
terminal.  This  problem  cannot  be  readily  solved  unless 
each  trace  is  recorded  to  disc  without  user  intervention. 
A*  p  r^s^nt ,  this  is  not  done  because  only  eight  traces  'two 
h-  :rs  f  data:  can  be  written  to  a  floppy  disc  before  it  is 


TIME  C  LOCRL ) 


TIME  RVEPfiGED  CT-2 


METERS  RRNGE  IN  METERS 


K3*-- 


1 3  i  2  a 

TIME  '  LCCOL  ) 


TIME  RVEPmGED  CT 


NPS  POOF 


t;*E  *=ivEcM'3ED  rr~2 


UPS  POOP 


TIME  (LOCRL) 


TIME  AVERAGED  CT-2 


NFS  POOF 


50 

TIME  ( LOCRL 
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